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Q\ • Abstract 

By means of solving the Bethe-Salpeter equation with the minimal perturbative kernel 
(ladder approximation) in the arbitrary covariant gauge, there was found a series of inter- 
O ■ cepts of meson Regge trajectories associated with colorless singularities of the four-quark 

Green function in SU(iV c ) quantum chromodynamics with massless quarks. The gauge 
and renorm-invariance of the result is grounded. 
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Introduction 



Q\ '. One of the most important and interesting problems of the theory of strong interaction is the 

problem of calculation of Regge trajectories (singularities of scattering amplitudes continued 
analytically to the region of complex angular momentum [TJ) in the framework of quantum 
^ ■ chromodynamics (QCD). Upon having the analytical behavior of Regge trajectories determined 
we automatically get information on hadron masses and decay widths corresponding to those 
points of Regge trajectories where trajectories themselves take on integer non-negative values. 
Notably, mass M and decay width T of some spin-iV hadronic state placed on Regge trajectory 
a{t) obey the equation 

a(t N ) = N (Ret N = M 2 , lmt N = -MT). 

Moreover, the behavior of leading Regge trajectories in the region of small negative values of 
the argument governs the high energy evolution of the diffractive pattern at scattering to small 
angles. 

From the very nature of Regge trajectories in quantum field theory (poles of scattering am- 
plitudes and corresponding Green functions) there follows one of their most important prop- 
erties — renorm-invariance, i.e., their independence of the renormalization scheme and the 
renormalization scale. 

The most popular approach to calculation of Regge trajectories in QCD is the Balitsky- 
Fadin-Kuraev-Lipatov (BFKL) approach, which, in fact, brings us to looking for poles of so- 
lution of the linear integral BFKL equation [2] - the modification of the Bethe-Salpeter (BS) 
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equation where external particles (quarks or gluons) and their propagators are replaced with 
the so-called "reggeized" partons, i.e., reggeons with quantum numbers of quarks and gluons 
obtained via solving the BS-equation with some perturbative kernel. Under usage of the BFKL 
approach the asymptotic behavior of some Regge trajectories was determined in the region 
of high negative values of the argument [31 H] where perturbative techniques are applicable. 
However, the problem of calculating renorm-invariant Regge trajectories in the region of small 
values of the argument (or, at least, renorm-invariant intercepts of Regge trajectories) is so 
far unsolved. The calculation of the "hard pomeron" (the leading reggeon associated with the 
spectrum of states formed by two "reggeized" gluons) intercept in the framework of the BFKL 
approach brings to the result explicitly depending on both the renormalization scheme and the 
renormalization scale [5]: 



where fj, is the renormalization scale and a s (fi) = S^i^- is the QCD effective running coupling. 

The direct aim of this paper is the calculation of renorm-invariant intercepts of Regge 
trajectories associated with colorless hadronic states formed by quark and antiquark. For this 
purpose we will use a technique alternative to the BFKL method that was applied earlier by 
C. Lovelace to the scalar field model 0g [6]. This technique consists in looking for singularities 
of the solution of the BS-equation with some perturbative kernel and implies usage of effective 
running coupling. The basic physical idea follows. The Green function singularity points 
associated with the complex value of the quark- antiquark system invariant mass should not 
depend on the quark virtualities (this is the main hypothesis based on the fact that Regge 
trajectories are functions of one dynamical variable). So we may give virtualities of the quark 
and the antiquark asymptotically high negative values (but the value of the system invariant 
mass itself can be put low or even equal to zero). Considering the BS-equation in the region of 
asymptotically high quark virtualities allows us to choose the renormalization scale large enough 
for having a possibility to use the approximation of perturbative kernel minimal on powers of 
the running coupling (instead of the full renorm-invariant kernel). Such an approximation 
is called the "ladder" one. At this, the running coupling itself demonstrates characteristic 
perturbative behavior at the chosen renormalization scale (connected with high virtualities of 
the quarks). The positions of Green function singularities must depend on the quark-antiquark 
system invariant mass and QCD fundamental parameters only. Renorm-invariance should reveal 
itself in the fact that calculated Regge trajectories must be independent of the renormalization 
scheme. 

The ladder approximation and intercepts of meson Regge 
trajectories 

In the case of zero invariant mass, t = 0, the BS-equation for the full four-quark Green function 
G(q,p; /i;g s (n)) (here q and p are 4-momentums of quarks) cut at two points is of the form 
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For the full renorm-invariant kernel K(q,p; [i\ <7 S (/•*)) there takes place a formal expansion 
on powers of the effective running coupling: 
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If there is a hard dynamical scale in the BS-equation (for example, large virtuality of quarks), 
then by proper choice of the renormalization scale one can make terms of order O(pf) in the 
kernel expansion to be little enough (in comparison with the tree term) to have a possibility to 
neglect them. We are interested in those singularities which depend on the invariant mass t of 
the system "quark-antiquark" and do not depend on other dynamical scales. So for considera- 
tion we may choose the region of large quark virtualities. At this, there emerges a possibility 
of usage of the minimal perturbative kernel and perturbative expressions for effective running 
parameters at scales of order of the quark virtualities. 

In such an approximation the connected part of G turns out to be equal to the sum of the 
ladder diagrams: 




Since we deal with the sum of the ladder diagrams (instead of the full Green function) 
the gauge and renorm-invariance of the results are not a priori (in contrast to the full kernel, 
the minimal perturbative kernel is not gauge and renorm-invariant) and should be grounded 
additionally. We will start from the arbitrary covariant gauge. 

In this paper we restrict ourselves by looking for the colorless singularities at zero invariant 
mass, t = (intercepts of meson Regge trajectories), in the case of massless quarks. Upon 
the average over quark colors in the colorless t-channel, the BS-equation for G in the ladder 
approximation takes the form of (for better visual perception we omit spinor indices at matrix 
structures) 

G{q,p)=i{2^5\q-p) + 

where q, p, and k are quark 4-momentums, 7° are Dirac matrices, q = g M 7 M , N c is the number 
of colors, g s {n) is the QCD effective running coupling, is the running gauge parameter, 
and n = \J—q 2 is the chosen renormalization scale (here and in the following, matrix structures 
to the left of the Green function should be attributed to the lower spinor line in the figures, 
and matrix structures to the right of G - to the upper one). 

Note that for large enough values of the renormalization scale the running gauge parameter 
behaves as £(//) ~ —5- ~ g 2 s {^) 0. So since we neglect loop corrections of order and higher 

in the kernel, we may also neglect the longitudinal part of the gluon propagator in the tree 
term of the kernel because it also gives a contribution of order g^. In other words, in the region 



of asymptotically high ^/—q 2 the approximation of the minimal perturbative kernel (ladder 
approximation) does not contain dependence on the gauge. 

After the convolution multiplication by q x q (the sign "x" denotes the direct product) we 
come to the equation 

qaqAl a G{q,p)l P ] = *(27r)V(g - p)q a qp[j a x /] + 
iV c 2 -l^ 2 (V^) f d*k ( (q-k) a (q-k) p \ a * & 

+ 2N C i(2*y J J^W \ 9a " (i-ky ) b G ( k > p W- a) 

Any ladder diagram can be represented as the sum of terms containing the direct product of 
two matrix structures. Since we put quarks massless, then each of these structures corresponds 
to the convolution product of the even number of Dirac matrices (vertices and quark propagators 
supply individual Dirac matrices to the corresponding structure). So G(q, p) may be represented 
in form of 

G(q,p) = G 00 (q,p)[I x /] + l -(G w (q,p)(q,p x - p,q x ) + E 20 (q, p)e^\s ri q S p n ) [a^ x /] + 

+ 1 -(G 02 (q,p)(q u p e - Pu q e ) + F 02 (q,p)e^ T q^)[I X a^} + 

+ ^(G 22 (q,p)(q^Px - Pnq\){quPe - Puqe) + E 22 (q,p)e^ X s v q 5 p n (quPe - Puqe) + 
+ F 22 (q,p)(q^p x - P^qx)e ue ^p T + H 22 (q,p)e^6 v q 5 p v e v ^p T )[a^ x a ve \ + (2) 
+ \(G 2 \q,p)(q, P x - p,qx) + E 2 \q,p)e flX s v q 5 p r, )[^ X X 7 5 ] + 

+l(G 43 (q,p)(q vPe -p v q e ) + F 42 (g, p)e^p T ) [ 7 5 X a^} + 

^40 „\L,5 ^ rl i ^04 / „ur xx Kl i r<44 



+G 4U (g,p)[ 7 5 x I] + G»\q,p)[I x 7 5 ] + G 44 (g,p)[ 7 5 x 7 5 ]- 



In square brackets there are Dirac matrix structures (J is the unit matrix, <j pv = i 1 ^ 1 ~ 7 7<J , 
7 5 = i7°7 1 7 2 7 3 ) and all unknown functions are scalars. 

After the substitution of the last expression into (PQ) we convolve it with the structure 
^[7 P x 7 " + D p x D a \ (here D p = ry^ 5 ) over all spinor indices. 
In view of relations 

q«qp Sp\Yl a G{q, p) 7 V + D^<*G(q, p)^D°] = q 2 G(q, p) 

lb 

and 

ff - (g ~^~ fe) ^ W7 a G(?, p) 7 V + D p ^ a G(k, ptflT) = 3 G(k,p) , 

where G(q,p) = G 00 (q,p) + (G 22 (q,p) - H 22 (q,p))(p 2 q 2 - (pq) 2 ) + G 44 (g,p), we come to the 
scalar integral equation 

e «(,, P ) = w *< 9 - P ) 5 > + / (^(m (3) 

for some linear combination of structure functions, G(q,p). In other words, we have imple- 
mented a partial diagonalization of Eq. (pQ). 



If we put g1{-\/—q 2 ) = 77i 7^-7 — (where n/ is the full number of quark flavors), 

Eq. ([3]) can be solved directly. In this way we can obtain some part of the spectrum of 
intercepts of meson Regge trajectories. The one-loop approximation of the running coupling 
cjsin) is not renorm-invariant. It contains the QCD dimensional parameter A depending on the 
renormalization scheme. Also, its value explicitly depends on the renormalization scale (above 
we have chosen /1 = \/—q 2 )- But we are interested in those singularities of G(q,p) which are 
renorm-invariant, i.e. independent of the renormalization scale (quark virtualities) and A. Only 
such singularities should be considered as intercepts of Regge trajectories. 

The procedure of solving ((3]) is expounded in the Appendix. Here we give only the obtained 
series of intercepts 

a W (0 ) = 9 ^ ~ V 1 (4) 

where k is an arbitrary integer and n/ is the full number of quark flavors. 



Discussion 

Let us point out the main features of series (HI). The calculated intercepts depend on the 
fundamental dimensionless parameters, N c and rif, and do not depend on quark virtualities, 
dimensional QCD parameter A, and the coupling. They do not depend on the renormalization 
scheme either. 

The spectrum of intercepts from (Jl]) will not change if we multiply the renormalization 
scale /1 (fixed in the BS-equation) by arbitrary constant, i.e., if we put /i = C^f—q 2 instead 
of fi = \f—q 2 - This is the straight consequence of their independence of q 2 , p 2 , and A. It 
is curious that independence of the coupling directly follows from the renorm-invariance of 
Regge trajectories in asymptotically free massless field models [8]. All these facts point to 
renorm-invariance of series (Jl]). The gauge invariance follows from the independence of the 
minimal perturbative kernel of the gauge parameter (this was discussed above). We would 
like to emphasize that the obtained result (for the case of massless quarks) is asymptotically 
accurate since the value of — q 2 is chosen to be asymptotically high. 

In our calculations the masslessness of quarks (chiral symmetry) was an essential point. In 
reality the chiral symmetry is broken, either explicitly (in the Lagrangian) or dynamically, and 
the question how spectrum (jl]) is affected by nonzero quark masses is very important. Mass 
terms are present in both the quark propagators and the effective running coupling. However, 
even the heaviest quark mass does not bias the asymptotical behavior of the effective running 
coupling in the region of extra-high values of the renormalization scale. Relative to quark 
propagators the situation is more complicated. For massive quarks the BS-equation takes the 
form 

(q - mi)G{q,p){q - m 2 ) = i(2vr) 4 <5 4 (g - p)[(q - mi) x (q - m 2 )\ + A(q,p) , 

where mi and mi are the quark masses and A(q,p) is the second (integral) term in the right- 
hand side of Eq. ([T]). In this case the numbers of both the Dirac matrix structures and the 
independent structure functions in (J2]) appreciably increase and partial diagonalization of the 
BS-equation cannot be implemented in such a simple way as we do for deriving ([3]) from ([I]). So 
at the moment we are not able to give quantitative estimation of the chiral symmetry breaking 
influence on series (jl]) explicitly from the BS-equation. But it is clear that increasing of quark 
masses should increase masses of meson states placed on Regge trajectories and, hence, decrease 
corresponding intercepts. Below we will see that comparison of (jl]) with meson phenomenology 
allows us to estimate this decrease for the leading {k = 0) intercept from (jl]). Besides, such 



decreasing of intercepts implies that spectrum (jj]) can be associated only with those observed 
reggeon families which have all intercepts lower than the leading intercept from (pi}. 

Let us consider light (i.e., composed of u- and d-quarks) mesons with isospin 1 = 1. This 
choice is motivated by the wish to avoid mixing with more heavy quark flavors. For such 
mesons there exist four families of Regge trajectories - a-reggeons, p-reggeons, 7r-reggeons, and 
6-reggeons. Chew-Frautschi plots (dependences of resonance spins on their masses squared) 
for these reggeons demonstrate unified linear behavioi0 with slopes about 1 GeV -2 [1]. At 
N c = 3 and rif = 6 (three generations of fermions in the Standard Model) the leading intercept 
from ([4]) is afq(0) = 1/7. Even-spin state reggeons with positive intercepts should not have a 
spin-0 state because it would have imaginary mass (tachyon). The leading 7r-trajectory has a 
spin-0 state, pion. Hence, if we assume the existence of the continious chiral limit, spectrum 
(jlj) cannot be associated with 7r-reggeons. The leading intercepts of p- and a-reggeons can be 
estimated not only via linear continuation of their Chew-Frautschi plots but also from the high 
energy evolution of angular distributions for reactions ir~ + p — > ir° + n and 7r _ + p — > r] + n. 
They both have values of about 0.45 [9|. Consequently, series @ can be associated only with 
b- reggeons. In Table [TJ we compare intercepts from series flU) at N c = 3 and rif = 6 with the 
observed states of 6-reggeons. 
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Table 1: The comparison of intercepts from (j4j) for the case N c = 3 and ri/ = 6 with the observed 
states of 6-reggeons. The subscripts denote resonance spins and the numbers in parentheses 
denote approximate values of meson masses in MeV. 
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Figure 1: Expected behavior of odd-spin meson Regge trajectories with negative intercepts. 

Series (jl]) accumulates towards the value ot^ifS) = —1 (similar accumulations take place 
for the four-gluon Green function [TQj ) . However, the intercept of the Chew-Frautschi plot 
corresponding to the next to leading (k = 1) 6-trajectory is lower than —1. The same takes 

1 This, of course, does not imply that true Regge trajectories are strictly linear since equations like a(M 2 — 
iAiT) = N do not imply that Rea(M 2 ) = N and Ima(M 2 ) = 0. Unfortunately, for most of the reggeons the 
only way to estimate their intercepts phenomenologically is the continuation of Chew-Frautschi plots. 



place for p-reggeons: Chew-Frautschi plots for all observed p-trajectories, except the leading 
one, have intercepts lower than —1 (the list of all observed p- and other mesons can be found at 
[TT]). But any true odd-spin Regge trajectory can take value —1 nowhere because this would 
cause a scattering amplitude unphysical singularity associated with some negative-spin state. 
Consequently, if we assume monotony of Regge trajectories below the first threshold, then all 
odd-spin Regge trajectories with negative intercepts should behave as in Fig{TJ Such a pattern 
seems quite natural. The nonlinearity of trajectories gains strength with radial quantum number 
k increasing, and true values of reggeon intercepts should have values higher than those of their 
Chew-Frautschi plots. 

Now we can estimate the degree of influence of chiral symmetry breaking on the leading 
intercept value from series (J4]). The true intercept of the leading b- reggeon must have value 
lower than afg(0) = 1/7 and higher than the intercept of the corresponding Chew-Frautschi 
plot, ajjn = —0.17 ± 0.04. For other 6-reggeons we cannot do such a concrete estimation, but 
it is clear that they should have intercept values higher than —1 and the accumulation of these 
intercepts as in is preserved at chiral symmetry breaking. So there takes place a qualitative 
correlation between the obtained series (j3J) and the phenomenology of 6-mesons. 

At the very end we must note that, in fact, we have not solved the BS-equation (pQ) com- 
pletely and found only that part of the spectrum of intercepts of meson Regge trajectories which 
corresponds to 6-reggeons. The full solution of Eq. (JT|) should also contain other parts of the 
four-quark Green function singularities which could be associated with a-, p-, and 7r-reggeons. 
Together with the problem of the chiral symmetry breaking influence on reggeon spectrums 
there is much further work to be done on these questions. 
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Appendix 

We will solve the integral equation 

e &,, P) = W v h -,W + ^^§^ / (^<5 (M (5) 

using the technique applied earlier to the BS-equation in the scalar field model [5]. 
After the Wick rotation Eq. (j3j) turns into 

mm - <*w?- *>?■ + 4^ w 1 / • (6) 

Similarly to the three-dimensional case the amplitude G can be expanded into a series of 
partial harmonics relative to the scattering angle 

a (PQ) 
cos v — 



\p\ l<?1 



A more detailed description of this formalism can be found in [6j [12] • We restrict ourselves to 
the list of basic formulas. 



In four space dimensions the partial wave expansion is of the form 

oo 

T(cos9) = Y / (2L + 2)U L (co8B)T L , 

L=0 

where Ul(x) are the second kind Chebyshev polynomials. The reverse formula is 

1 rl 

t l = —r T tt / dxVl - x 2 U L (x) T(x) . 

7T(L + 1) J-l 

If scalar functions A, B, and C are related by 

A(q,p) = J^yj dkB(q,k)C(k,p) } 
then for the partial wave projections (in what follows, q = \q\, k = \k\, p = 

I roo 

A L (q,p) = —j o dkk 3 B L (q,k)C L (k,p). 
And, at last, the definition of the four- dimensional analog of the second-kind Legendre function 



is 



Ql(z) = — — - / dxVi 



tt(L + 1) J-i z — x 

For this function the following integral representation exists: 

-(L+l)|r| 
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Q[ 4) (coshr) = J°°dT'e-^ 
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Now we introduce a dimensionless function G^Qnq, hip) = G . L 2 q f^ ■ Then, after proceeding 



in (El) to the partial wave representation we come to 

A?-lStf( 9 ) f^,,^^J 

We define the Fourier transform of G^flng, Inp) as 

I roo roo 



G L (ln q, hyp) = <?(ln g - lnp) + ^ 1 3 9 ^ Jd(lnk)Q^ (cosh (in-)) G L (ln k, hip) . (7) 



rOO rOO _ 

G L (lng,lnp) = — / db db'e i{blnq - b lnp) [(L + l) 2 + b 2 ]F L (b,b') 

ZTY J —oo J —oo 



where 



1 roo roo „ 

F L (\nq, \np) = — / d& / db'e i{blnq - b lnp) F L (b, b') . 

2,71 J-OO J-OO 

Taking into account that 

£^^(coshr) = (L+1 2 )2 + fc2 

we come to the equation 

1 roo roo sis ~ 

— / db db'e i{blnq - b lnp) [(L + l) 2 + b 2 }F L (b,b') = 5{lnq - lnp)+ 

27T ,7 — oo J —oo 

N 2 -13f s (q) l_ roo roo i(blnq _ b , lap)p 
2N C An 2 27ri-oo i-oo K J 



which is, in turn, equivalent to the differential equation 

d 2 



(L+l) 
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/V 2 — IS o 2 ^ 

F L (lng, hip) = 5(lng - lnp) + |^F L (lng, hip) . (8) 



In what follows, at asymptotically high values of q we will use the perturbative one-loop ap- 



proximation to the QCD effective running coupling, g 2 (q) 
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(at this, all the 



dependence on the renormalization scheme is in the constant A). So the Fourier transform of 
(IH1) takes the form of 
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[(L + l) 2 + b 2 }F L (b, b') = 5{b ~b')-— db"e(b - b")F L (b", b' 



(9) 



where e(b) = b/\b\ and E 



9(^-1) 



N c (llN c -2nf) ' 

Now it is convenient to introduce a new variable, 



X 



db" 



arctan 



L+l 
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/o (L + l) 2 + 6" 2 L + l ' /vv " v 2(L + 1) 

Function x(b) maps 6 = (— oo, oo) onto x = ( — ^) monotonically. So 

e(b - V) = e( X - X') , 6{b - b'){{L + l) 2 + b' 2 ) = 6( X - X ') , 

where X ' = xip') and in terms of the new function yi(x, x') = ((L+l) 2 +b 2 )F L (b,b')((L+l) 2 +b' 2 ) 
Eq. takes the form of 



vl(x,x') = Kx-x!) 



iE r°° 



dx"<X-x")Ux'\x') 



After substitution yL(x,x') 



9 ^q^ - we come to differential equation 



Kx - x) + l -^{m,x') +j(-n,x')} = e ~ lEx ^ {e lEx Jix,x)} ■ 



The solution of this equation is 

e lEx j(x,x') = e lE *' 
The value of C is not arbitrary. Since 



C+^eix-x'] 



j(±n,x') = e 



then adding ± forms of this equation, we find C = |tan(ivO). Hence, finally, the solution of 
(Q is of the form 

[(L + l) 2 + b 2 ]F L (b,b') = 5(b - b') - l -^ e iEMb,) ~ x{b)) [e(x - x') + itan(£fi)] 

from where we obtain a series of the amplitude singularities (intercepts of meson Regge trajec- 
tories) 

9(N 2 - 1) 



4 } (o) 



(2k + 1)JV C (11JV C - 2n f ) 



- 1 



where k is an arbitrary integer. 
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